Cauchy nights in confining potentials 



43 



O 

> 



I. MOTIVATION 



Piotr GarbaczewskQ 
Institute of Physics, University of Opole, 45-052 Opole, Poland 

We analyze confining mechanisms for Levy flights evolving under an influence of external po- 
tentials. Given a stationary probability density function (pdf), we address the reverse engineering 
problem: design a jump-type stochastic process whose target pdf (eventually asymptotic) equals 
the preselected one. To this end, dynamically distinct jump-type processes can be employed. We 
demonstrate that one "targeted stochasticity" scenario involves Langevin systems with a symmetric 
stable noise. Another derives from the Levy-Schrodinger semigroup dynamics (closely linked with 
topologically induced super-diffusions), which has no standard Langevin representation. For compu- 
tational and visualization purposes, the Cauchy driver is employed to exemplify our considerations. 
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, We consider a subclass of so-called Levy nights that is mathematically identifiable as symmetric stable stochastic 
processes, [J- [II]. These non-Gaussian processes are of the jump-type and, in contrast to familiar Gaussian diffusion- 
type processes, involve a number of obstacles. 

One of them is a shortage of explicit analytic solutions, e.g. explicit probability density functions (pdfs) and 
transition probability densities. Another comes from theoretically admitted existence of arbitrarily small jumps and 
<"*> |- fat tails of the pdf which typically preclude the existence of (in the least the second) moments. In the presence of 
confining potentials the resultant pdfs may admit higher moments, but merely a finite number of them may exist. 
Third, some care is needed in any computer-assisted analysis of Levy Sights, since imposing a lower and upper bound 
on the size of jumps, sets the problem within the ramifications of the central limit theorem which implies the standard 
C Gaussian limit, [a], albeit the simulated process is non-Gaussian by definition. 

In the present paper we set general confinement criterions for symmetric stable processes, with special emphasis on 
| the analytically tractable case of the Cauchy noise With regard to a specific response to external potentials, 

we pay special attention to two classes of jump-type processes: those related to the Langevin equation and those 
[• — ' induced by the Cauchy-Schrodinger semigroup dynamics (that involves a fractional analog of the generalized diffusion 
i equation). We leave for a separate study another interesting possibility, that is based on Bochner's concept of 
subordination, c.f. (l~7j . 

Levy flights in confining potentials with a standard (mostly for an additive noise) Langevin representation have 
£"-~. ■ received ample attention, and [HI, For Cauchy-Langevin processes, a manipulation with the forward drift 

and/or its (external conservative force) potential sets rules of the game, e.g. directly leads to stationary probability 
densities. They never have a Gibbs-Boltzmann form, characteristic for Gaussian diffusion processes, jl7ll . 

Another class of Levy processes, that are driven by dynamical semigroups, was analyzed in detail in Cauchy 



driver has received there special attention. The semigroup-driven processes independently reappea red (in the context 
of systems with topological complexity like folded polymers or complex networks) in Refs. [20t]— [23l] . In Refs. [TsL f]~9l] . 
external potentials appear as additive perturbations of the noise generator and under the name of effective potentials 
they appear in [2fj| as well. In the "topological" literature, the resultant semigroup dynamics has been implemented 
via local modifications of jump rates of the associated jump-type process. 

The semigroup-driven (and topologically-induced) processes appear to have no standard Langevin representation, 
irrespective of whether we adopt an additive or multiplicative noise. Only under special circumstances a connection 
with the multiplicative noise has been established in Rcf. [201 ] . but it is not a generic property of semigroup-driven 
processes, see e.g. also 25? 1. 

Although we formulate a framework incorporating symmetric stable processes in general, we strongly rely on a 
mathematical theory of the Cauchy semigroup-driven processes. This theory, without any "topological" context, has 
been formulated in Ref. No explicit analytic examples of confining potentials, nor pdfs were given there. 

We shall demonstrate that the Langevin-driven and semigroup-driven Levy processes stay in affinity and may share 
common for both stationary pdf. A super-diffusive dynamical pattern of behavior is generically expected to arise. An 
asymptotic approach towards a stationary pdf is then in principle possible. 
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This motivates our " targeted stochasticity" discussion whose original formulation for Langevin-driven Levy systems 
can be found in Ref. [l7|, the reverse engineering problem being included. The necessity of considering the semigroup 
dynamics in this reconstruction problem points first, to an independent dynamical mechanism, and second, to the 
fact that only the properly tailored semigroup (and thus the semigroup potential) may guarantee that the prescribed 
invariant density actually is an asymptotic target of the process. This point was left aside in "topological" references 
[13]- [23|], where the pertinent invariant densities were postulated to have a Gibbsian form, hence a tacit assumption 
was made about extyremely strong confining properties of the topologically-induced process. 

The original reverse engineering problem reads: given a stationary pdf, can we tailor a drift function so that the 
system Langevin dynamics would admit the predefined as an asymptotic target ? In the course of our discussion, 
we in fact extend the range of applicability of the original "targeted stochasticity" scenario and demonstrate that, 
for a symmetric stable driver, a priori chosen stationary pdf may serve as a target density for both Langevin and 
semigroup-driven jump-type processes. Even though their detailed dynamical patterns of behavior are different. In the 
near-equilibrium regime this dynamical distinction becomes immaterial. For analytically tractable and visualization 
insights, we shall basically refer to Cauchy processes in confining potentials. 

For the record we point out that the term "equilibrium" needs to be addressed with some care for non-Gaussian 
processes. No physical thermalization mechanisms have ever been proposed for Levy flights. Moreover, their physical 
"reason" (origin of noise) appears to be exterior to the physical system, with no reliable kinetic theory background, 
and therefore no fluctuation-dissipation response theory could have been been set for any stable noise. 

To the contrary, the noise "reason" is definitely an intrinsic feature of the environment-particle coupling in case 
of the standard Brownian motion, based on the kinetic theory derivations. All traditional fluctuation-dissipation 
relationships find their place in the Brownian framework. None of them has been reproduced in the context of Levy 
flights. 



II. RESPONSE OF GAUSSIAN NOISE TO CONFINING POTENTIALS: SMOLUCHOWSKI 

PROCESSES 

Albeit we are primarily interested in jump-type stochastic processes, certain useful intuitions can be borrowed from 
the standard theory of Brownian motion. Namely, let us consider a one-dimensional Smoluchowski diffusion process 
[2^1, with the Langevin representation 

x = b(x,t)+A(t) (f) 

where (A(s)) = 0, (A(s)A(s')) = 2DS(s — s') and b(x) is a forward drift of the process having the gradient form 
b = 2_DV<f>, where D stands for a diffusion constant. 

If an initial probability density po(x) is given, then the diffusion process obeys the Fokker-Planck equation 

dtp = DAp - V (b ■ p) . (2) 

We introduce an osmotic velocity field u = Dlnp, together with the current velocity field v = b — u. The latter 
obeys the continuity equation dtp — — Vj, where j — v ■ p has a standard interpretation of a probability current. 

Presently we pass to time-independent drifts of the diffusion process, that are induced by external (conservative, 
Newtonian) force fields / = — VV. One arrives at Smoluchowski diffusion processes by setting 

b = -^ = -^VV. (3) 

mp rap 

This expression accounts for the fully-fledged phase-space derivation of the spatial process, in the regime of large (3. 
It is taken for granted that the fluctuation-dissipation balance gives rise to the standard form D = kBT/m(3 of the 
diffusion coefficient. 

Let us consider a stationary asymptotic regime, where j — > = 0. We denote an (a priori assumed to exist) 
invariant density = p*(x). Since v„ — and b — f/m/3, by its very definition, does not depend functionally on the 
probability density, there holds 

b = 6* = = D V In p* . (4) 

Consequently, we have p*(x) = (\/Z) exp[— V(x)/fcsT], where 1/Z is a normalization constant. Our outcome has the 
familiar Gibbs-Boltzmann form. 

Following a standard procedure [26l . f20j ] we transform the Fokker-Planck equation into an associated Hermitian 
(Schrodinger-type) problem by means of a redefinition 

p(x,t)=V(x,t)pl /2 (x) (5) 
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that takes the Fokker-Plack equation into a parabolic one, often called a generalized diffusion equation: 

= L>A* - Vf . (6) 
Its potential V derives, as a function of the drift b(x), from a compatibility condition 

V(x) = (l/2)[6 2 /(2^) + V&]. (7) 

In view of Eq. (4), its equivalent form is 

V{x)=D^l. (8) 
p* 

If the (l/2mD rescaled) Schrodinger-type Hamiltonian H = —DA + V is a bounded from below, self-adjoint 
operator in a suitable Hilbert space, then one arrives at a dynamical semigroup exp(— tH), with the dynamical rule 
^>(x,t) — [exp(—tH)if?](x,0), pushing forward in time the initial data ^(x,0). The semigroup is contractive, hence 

1/2 

asymptotically ty(x,t) — + pj . Accordingly, p(x,t) — ► p*[x). 

The above Schrodinger semigroup (parabolic) reformulation of the Fokker-Planck equation refers to the very same 
diffusion process and the dynamics of p{x, t) does not depend on the theoretical framework of choice. In below we 
shall demonstrate that for non-Gaussian processes, the semigroup-driven and Langevin-induced dynamics refer to 
incquivalent dynamical patterns of behavior. Even if both are associated with a common stationary (target) pdf. 

III. RESPONSE OF LEVY FLIGHTS TO CONFINING POTENTIALS 

A. Levy driver 

Let us set general rules of the game with respect to the response of any symmetric stable noise to external potentials. 
We recall that a characteristic function of a random variable X completely determines a probability distribution of 
that variable. If this distribution admits a density p(x), we can write < exp(ipX) >= J R ,p(x) exp(ipx)dx which, for 
infinitely divisible probability laws, gives rise to the famous Levy-Khintchine formula (see, e.g. [3|]). 

From now on, we concentrate on the integral part of the Levy-Khintchine formula, which is responsible for arbitrary 
stochastic jump features: 



F(p) 



-Hoc. 



exp(ipy) - 1 - 

1 + y l 



"(dy), (9) 



where v(dy) stands for the appropriate Levy measure. The corresponding non-Gaussian Markov process is character- 
ized by < exp(ipX t ) >— exp[— tF(p)] and, upon setting p = — iV instead of p, yields an operator F(p) — H which is 
a direct analog of the free Schrodindger Hamiltonian. 

We restrict further considerations to non-Gaussian random variables whose probability densities are centered and 
symmetric, e.g. a subclass of stable distributions characterized by 

F(p) = \\p\^ ^> H = \\A\"/ 2 . (10) 

(In passing, we note that the adopted definition of a pseudo-differential operator may be replaced by the negative of 
a suitable Riesz fractional derivative.) In the above, p, < 2 and A > stands for the intensity parameter of the Levy 
process. The fractional Hamiltonian H , which is a non-local pseudo-differential operator, by construction is positive 
and self-adjoint on a properly tailored domain. A sufficient and necessary condition for both these properties to hold 
true is that the pdf of the Levy process is symmetric, [3j. 

The associated jump-type dynamics is interpreted in terms of Levy flights. In particular 

F(p) - AH - H = F(p) = A|V| = A(-A) 1 / 2 (11) 

refers to the Cauchy process, see e.g. |5, ITR [19} ■ 

The pseudo-differential Fokker-Planck equation, which corresponds to the fractional Hamiltonian (fTU)) and the 
fractional semigroup exp(— tH^) — exp(— A| A^/ 2 ), reads 

d tP = -\\Ar/ 2 p, (12) 
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to be compared with the Fokker-Planck equation for freely diffusing particle dtp = DAp. 

For a pseudo-differential operator |A| M ' 2 , the action on a function from its domain is greatly simplified, in view of 
the properties of the Levy measure v^dx). Namely, remembering that we overcome a singularity at by means of 
the principal value of the integral, we have [HI, [lj| : 

(\A\^f)(x) = - [ [f(x +y)- f(x)]^(dy) . (13) 



By changing an integration variable y to z = x + y and employing a direct connection with the Riesz fractional 
derivative of the p-th order, [T3 |. we arrive at 

(|Ar/2 f)(x) = 1)^(^/ 2 ) i m-m dz (14) 

with (|A^/ 2 /)(x) = -d fl f(x)/d\x\f. The case of fi = 1 refers to the Cauchy driver (e.g. noise). 

We note a systematic sign difference between our notation for a pseudo-differential operator |A|^/ 2 and this based 
on the fractional derivative notion, like e.g. A^/ 2 = d^/dlx^ of Refs. [20L |2^|. 



B. Langevin scenario 



In case of jump-type (Levy) processes a response to external perturbations by conservative force fields appears 
to be particularly interesting. On the one hand, one encounters a widely accepted reasoning (Refs. @-[I3]) where 
the Langevin equation, with additive deterministic and Levy "white noise" terms, is found to imply a fractional 
Fokker-Planck equation, whose form faithfully parallels the Brownian version, e.g. (c.f. [|, see also ||) 

x = b{x)+A> 1 (t) 



d tP =-\7(b-p)-X\A\^ 2 p. (15) 

We emphasize a difference in sign in the second term, if compared with Eq. (4) of Ref. [f|. There, the minus sign is 
absorbed in the adopted definition of the (Riesz) fractional derivative. Apart from the formal resemblance of operator 
symbols, we do not directly employ fractional derivatives in our formalism. 

Let us assume that the fractional Fokker-Planck equation (|15[) admits a stationary pdf p*(x). Then, a functional 
form of the drift b(x) can be reconstructed by means of an indefinite integral 

J\A\^p4x)dx 

b(x) = -A- — . (16) 

p*{x) 

This is the reverse engineering problem of [T3|. 



C. Levy-Schrodinger semigroup 

On the other hand, by mimicking the previous Gaussian strategy, we can directly refer to the Hamiltonian framework 
and dynamical semigroups with Levy generators being additively perturbed by a suitable potential, see e.g. [il^.[l9l|. 
For example, assuming that the functional form of V(x) guarantees that = A| A| M//2 + V is self-adjoint and bounded 
from below in a suitable Hilbert space, we may readily pass to the fractional (non-Gaussian, jump process) analog of 
the generalized diffusion equation: 

d t ^ = -A|A|' l/2 * - V*. (17) 

The dynamical semigroup reads exp(— tH^) and the compatibility condition affine to that of Eq. ([8]), typically takes 
the form of the time-adjoint equation for an auxiliary function 9(x, t): 



d t 9 = \\A\» /2 d + V9. 



(18) 
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General theory [HGJ] tells us that a (properly normalized) product 9*(x, t)6(x, t) determines a probability density 
p(x, t) of a Markov process that interpolates between the boundary data p(x, 0) and p(x, T), in the time span t 6 [0, T] . 
Let us assume that there exists an invariant (stationary) pdf p* of this Cauchy semigroup-induced process. If 

1 /2 

we demand that 8(x,t) actually does not depend on time, and adopt a decomposition p = \&p* , c.f. Eq. (5), we 

1/2 

are allowed to set 6 = p* and remove limitations upon the time interval. Ultimately, we arrive at a compatibility 
condition that is a direct fractional version of Eq. (J5J): 

V = -A ' ' J* • (19) 
p* 

1 /2 

In. the present case, we can readily evaluate the dynamics of p{x,t) — ^(x,t)pj (x): 

dtp = pl /2 d t ^ = ~\pl /2 \A\^ 2 [ P : 1/2 P } + v- P . (20) 

This is a departure point for the reverse engineering procedure: given the stationary pdf, find the semigroup potential 
V(x), 

1/2 

It is interesting to observe that by making cosmetic changes: set A = 1, formally identify pj = exp[— (3V(x)], with 
whatever V and (3 = 1/ksT, we and up with a familiar form of the transport equation previously introduced in a 
number of papers: 

dtp = - exp(-/?V/2) |A|^ /2 exp(/?U/2)p + pexp(/3U/2)|A|^ /2 exp(-/3U/2) , (21) 

c.f. formula (6) in [22| . formula (5) in [23| and formula (36) in [2(j. There, the investigated process was named a topo- 
logically induced super-diffusion. Wepoint out a systematic sign difference between our | A| M / 2 and the corresponding 
fractional derivative A' 1 / 2 of 0, [22|, [23| . Graphically these symbols look similar, but have different origin. 

Remark: It is of some interest to invoke an independent approach of Refs. [2l|,[22| where one modifies jumping rates 
by suitable local factors, to arrive at a response mechanism that is characteristic of the previously outlined semigroup 
dynamics. In view of (14), the free transport equation dtp = — XlA^^p can be re- written as a master equation 
dtp{x) — J[w(x\z)p(z) — w(z\x)p(x)}v f j_(dz). The jump rate w(x\y) ~ l/|x— is an even function, w(x\z) —w{x\z). 
If we replace the jump rate w(x\y) of the free fractional dynamics by the expression w^{x\y) ~ cxp[$^)-<^(y)] 

account for the fact that w<t,{x\z) ^ w c i > (z\x), then the master equation takes the form: (\/\)d t p = \/S.\^ 2 f = 
— (exp $) | A|^/ 2 [exp(— $)p] +pexp(— $)| A|^/ 2 exp(<I>). Whatever <I>(a;) has been chosen (up to a normalization factor), 
then formally = exp 2$ is a stationary solution of that transport eqution. We note that a physically attractive 
point in the topologically-induccd dynamics pattern was an assumption that exp 2$ sets a Gibbsian form of the pdf. 
Accounting for the normalization factor 1/Z one presumes that p» = (1 / Z)exp(—V* / hsT) with an external potential 



V„ = — fcsTln(Z p»), whose physical origin is based on a crude phenomenology, c.f. [2l|,|22j. With these re-definitions, 
the above transport equation takes the form (21). 

IV. REVERSE ENGINEERING FOR CAUCHY FLIGHTS 

By choosing p = 1 m the above, we narrow down the whole discussion to Cauchy processes, when e.g. (all integrals 
are evaluated by means of their Cauchy principal value) 

(lAI^/Xz) = - / [f( x + y)-f( x )-y^I^]u ,(dy) 

Jr. i + y 2 



(\Ar /2 f)(x) = / [f(x + y)- f(x)]^(dy) . (22) 



The Cauchy-Levy measure, associated with the Cauchy semigroup generator IAI 1 / 2 = |V|, reads 

vx/2{dy) = -%. (23) 
7T y z 
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By changing an integration variable y 



x + y, we give Eq. (20) the familiar form 



(|V|/)W = ~nj \z-x\ 2 dz 
where l/ir\z — x\ 2 has an interpretation of an intensity with which jumps of the size \z — x\ occur. 



(24) 



A. Ornstein-Uhlenbeck process 

In case of the Ornstein-Uhlenbeck-Cauchy (OUC) process, the drift is given by b(x) 
invariant density associated with 



reads: 



d tP = -\\v\p + v[{ 1 x)p] 



p*{x) 



n(a 2 + x 2 ) 



-jx, and an asymptotic 
(25) 

(26) 



where a = A/7, c.f. Eq. (9) in Ref. @ 

A characteristic function of this density reads —F(p) = — a\p\ and gives account of a non-thermal fluctuation- 
dissipation balance. The modified noise intensity parameter tr is a ratio of an intensity parameter A of the free 
Cauchy noise and of the friction coefficient 7. 

From the start we know what is the drift b{x) — —\x which directs the process towards a target (stationary) pdf. 
To deduce the Feynman-Kac potential V for the OUC process, we need to evaluate 



A 



2\l/2 



V{x) = -{a 2 +x 2 ) 



In the notation a = a 2 + x 2 , b = 2x, R(y) = <r 2 



a/ct 2 + (x + y) 2 y/a 2 + x 2 



dy 



(27) 



+ y) 2 indefinite integral reads, [3~i| 



y 2 ^W) V 2 



dy 



yy/a 2a \ 2a\y\ 



(28) 



Because of the singularity at y — 0, we must handle the integral in terms of its principal value, i.e. by resorting to 



J — » J_ + J e , and next performing the e — > limit. 

Taking into account that arsinhx = ln(x + y/1 + x 2 ), [13], we ultimately get 



V(x) 



■ In ■ 



la + x 



(29) 



V(x) is bounded both from below and above, with the asymptotics (2/|a;|) In \x\ at infinities, well fitting to the general 
mathematical construction of (topological) Cauchy processes in external potentials, [l9j . 

Accordingly, we know for sure that there exists a jump-type process driven by the Cauchy semigroup with the 
potential function V, Eq. (29), whose invariant density coincides with that for the Langevin-supported OUC process. 
This form of the semigroup potential, gives a guarantee that p* actually is an asymptotic invariant density of the 
process. In Fig. 1 we reproduce the functional shape of the potential (29), [33| . 



B. Confined Cauchy process 



The OUC process is not confined, since for the Cauchy density its second moment is nonexistent. We shall adopt 
the OUC discussion to Cauchy-type processes whose invariant densities admit the second moment. Let us consider 
the quadratic Cauchy density: 



p*{x) 



2 1 

7T (1 + z 2 ) 2 



(30) 
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FIG. 1: Upper panel: the coordinate dependence of semigroup potentials V(x): (29) for different a and (33) (inset). Lower 
panel: time dependent variance X 2 (t) for Lan gevin-type (solid line) and semigroup-driven (dashed line) processes associated 
with the pre-defined target pdf (33). Points correspond to numerical calculation. 



The action of |V| upon this density can be evaluated by recourse to the free Cauchy evolution. 

We note that (\/ \f2Tt)p 1 J 2 = (l/ir)/(l + x 2 ) actually is the Cauchy probability density. Let us consider f(x) = (lJ 2 
as the initial data for the free Cauchy evolution d t f = A|V|/. This takes f(x) into 

nx ' t) = n[(i + xtr+^}- (31) 



Since 



A|V|/ = -Hma t / (32) 
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we end up with 

w \ limtxo^t/ x 2 -1 

V{x) = (x) = A . (33) 

/ x z + 1 

A minimum — A is achieved at x = 0, V = occurs for x = ±1, a maximum +A is reached at x — > ±oo. The functional 
shape of the potential is depicted as an inset in Fig. 1. 

The potential is bounded both from below and above, hence can trivially be made non- negative (add A). Therefore, 
the invariant density (30) is fully compatible with the general construction of the Cauchy-Schrodinger semigroup and 
the induced jump-type process, c.f. Corollary 2, pp. 1071 in |l9j. This topological Cauchy process is induced by the 
Cauchy generator plus a potential function V given by Eq. (49). c.f. Corollary 2, pp. 1071 in [l!|. The process is of 
the jump-type and can be obtained as an e { limit of a step process, e.g. jump process whose jump size is bounded 
from below by e > but unbounded from above. 

In connection with the reverse engineering problem of Ref. [l7| let us note that if the quadratic Cauchy density 
actually stands for a stationary density of the fractional Fokkcr-Planck equation with a drift, then we should have: 

d t p* = 0= -V(6p*)-7|V|p.. (34) 
Therefore the drift function may be deduced by means of an indefinite integral: 

b{x) = [ (|V|p*)(x) dx = ~ \ (x 3 + 3x) . (35) 

For clarity of discussion, in the lower panel of Fig. 1, we report a comparison of dynamical patterns of behav- 
ior for the semigroup-driven and Langevin-induced scenarios beginning from common (delta-type) initial data and 
approaching a common (pre-defined) target pdf (33). 

Remark: It is worthwhile to notice that potential functions for the drift b = — V£7 of the form U(x) = ax 2 + bx A 
were investigated in Ref. [HI ] with a focus on bimodality of the resultant stationary pdfs. For example a quartic 
potential U — bx A is known to induce p sr {x) = [ir(l — x 2 + .t 4 )] -1 . Our unimodal density (33) belongs to the subclass 
associated with the just mentioned U{x). 



C. Confined Cauchy family 

We may consider various probability densities as trial ones. Let us pay attention to a broader class of densities that 
bear close affinity with the Cauchy noise. With a given continuous probability distribution p we associate its Shannon 
entropy S(p) = — J p \npdx, [32j. If an expectation value B[ln(l + x 2 )] is prescribed (e.g. fixed), the maximum 
entropy probability function belongs to a one-parameter family 

r(cy) i 

P * {X) = ^T{a- 1/2)) (l + x 2 )« (36) 

where a > 1/2, [11]. 

The gamma function r(a) = J Q exp(— t) t a ~ x dt we specialize to integer a = n + 1-values, with n > 0. Then 
r(n + 1) = nl and T(a - 1/2) ^ T(n + 1/2) = [{2n)\^\/n\2 2n . 
As an exemplary case let us consider 

p*(x) = -^- T / (37) 

57T (1 + X A Y 

By adopting the previous procedure, c.f. (3lj . and evaluating the principal value integrals, we end up with the following 
expression for the Cauchy semigroup potential: 

V(x)= 2(1 l x2) (x 4 + 6x 2 -3). (38) 

The potential is bounded from below, its minimum at x = equals —37/2. For large values of the potential 
behaves as ~ (j/2)x 2 i.e. shows up a harmonic behavior. 

Apart from the unbounded-ness of V(x) from above, this potential obeys the minimal requirements of Corollary 2 
in Ref. [l9j]: can be made positive (add a suitable constant), is locally bounded (e.g. is bounded on each compact 
set) and is measurable (e.g. can be arbitrarily well approximated by means of sequences of step functions). The 
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Cauchy generator plus the potential (38) determine uniquely an associated Markov process of the jump-type and its 
step process approximations. 

We can readily address the reverse engineering problem of Ref. [13] ■ For the density (37), we ultimately get, [33|: 

b(x) = - — (5x 6 + 2lx 4 + 35x 2 + 35) . (39) 
16 

This a bit discouraging expression which shows a linear friction b ~ —x for small x and a strong taming behavior 
b ~ —x 7 for large x, still fits to the above mentioned Corollary 2 of Ref. [l9j . 



V. SUMMARY 

We have generalized the reverse engineering (targeted stochasticity) problem of Ref. [TtJ beyond the original Levy- 
Langevin processes setting. We have demonstrated that the notion of Levy flights in confining potentials is not limited 
to the Langevin scenario. The Levy-Schrodinger semigroup involves the notion of external potentials as well. But 
then with no link to any standard Langevin representation. 

Our version of the reverse engineering problem amounts to reconstructing from a given (target) stationary density 
the potential functions that either: (i) define the forward drift of the Langevin process, or (ii) enter the Schrddinger- 
type Hamiltonian expression in the semigroup dynamics. Both dynamical scenarios are expected to yield the same 
asymptotic outcome i.e. the pre-selected target pdf. This goal can be achieved in the semigroup picture (and models of 
an impact of inhomogeneous environments upon Levy flights) only if suitable restrictions on the semigroup potentials 
are observed. The relevant mathemtical hints come from Ref. 19] and were illustrated for the case of Cauchy driver. 

Insightful, explicitly solvable models are scarce in theoretical studies of Levy flights, in the presence of external 
potentials and/or external conservative forces. Therefore, our major task was to find novel analytically tractable 
examples, that would shed some light on apparent discrepancies between dynamical patterns of behavior associated 
with two different fractional transport equations (15) and (20) that are met in the literature on Levy flights. Albeit 
the predominant part of this research is devoted to the standard Langevin modeling. 

We note that a departure point for our investigation was a familiar transformation of the Fokker-Planck operator 
into its Hermitian (Schrodinger-type) partner, undoubtedly valid in the Gaussian setting. The Fokker-Planck and the 
corresponding parabolic equation (plus a compatibility condition (5)) essentially describe the same random dynamics. 

An analogous transformation is non-existent for non-Gaussian processes. The two fractional Fokker-Planck equa- 
tions (13) and (15), are inequivalent in the non-Gaussian setting, hence the semigroup dynamics and the Langevin 
dynamics with the Levy driver (e.g. noise) refer to different random processes. This behavior we have depicted in 
the lower panel of Fig. 1. The main technical reason of the incongruence of the two processes seems to be rooted in 
that the stable noise generator is a (non-locally defined) pseudo-differential operator, while the standard Laplacian 
(Wiener noise generator) is locally defined. The reverse engineering problem allowed us to demonstrate that those 
two processes may nevertheless share the same target pdf and may interpolate between common pairs of boundary 
(initial and terminal) pdfs. Albeit in a different dynamical fashion. 

One may wonder whether there is some symmetry principle (like e.g. a remnant of the time-symmetric formulation 
of the Schrodinger boundary data problem, [3, E, H3, 13 that allows to relate two fixed boundary densities by 
means of different dynamical scenarios. Actually, our observation that the semigroup-driven and Langevin-driven 
jump-type processes may share a common invariant pdf, that in turn is dynamically accessible from a commmon 
for both processes initial pdf, stands for an indirect proof that the involved dynamical scenarios are different. The 
assumption about driving mechanisms is the only freedom left in the above mentioned boundary data problem, once 
the initial and terminal pdf data are chosen. 
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